Semiquantal dynamics of fluctuations: 
ostensible quantum chaos 
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The time-dependent variational principle using generalized Gaussian trial 
functions yields a finite dimensional approximation to the full quantum dy- 
namics and is used in many disciplines. It is shown how these 'semi-quantum' 
dynamics may be derived via the Ehrenfest theorem and recast as an ex- 

^ . tended classical gradient system with the fluctuation variables coupled to the 

O . 

■"sj" , average variables. An extended potential is constructed for a one- dimensional 

system. The semiquantal behavior is shown to be chaotic even though the 

'-q ■ system has regular classical behavior and the quantum behavior had been 

Q ■ assumed regular. 
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PACS number(s): 05.45. ,03.65.S,05.40. 



There has been substantial effort made, over the years, to understand quantum systems 



using a system of few classical variables. These can be motivated in the h — ► limit, such as 
effective potential techniques[l,2] and semi-classical methods like that of WKB, which lead 
to the Einstein-Brillouin-Keller (EBK) [3] quantization rules for integrable systems. These 
also include approximations to the Feynman path- integral formulation [4], used to derive 
the Periodic-Orbit Trace Formula for chaotic systems [5]. This relates the spectrum of the 
quantum system to a weighted sum over the unstable periodic orbits of the classical system. 
They can also arise, as in Quantum Chromo Dynamics for example [6], in the limit of large 
N (number of degrees of freedom), and a classical phase space can be shown to exist in the 
N = oo limit [7]. Further, there are many equivalent mean- field theories that are used in 
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nuclear physics, quantum chemistry, quantum field theory [8, 9, 10], condensed matter, statis- 
tial mechanics and optics (see the excellent reviews [11,12]). These are known variously as 
the time-dependent Hartree-Fock method, the Gaussian variational approximation, etc. One 
way to understand these approaches is to consider the time-dependent variational principle 
(TDVP) formulation [8,13,14] wherein one posits the action 



The requirement that ST = against independent variations of < and \&,t > yields the 
Schrodinger equation and its complex conjugate as the respective Euler-Lagrange equations 
(note that a ray rotation, i.e. \^,t >— > exp(iX(t)/h)\^,t > leaves the variational equations 
unchanged). The true solution may be approximated by restricting the choice of states to 
a subspace of the full Hilbert space and finding the path along which ST = within this 
subspace. In exceptional circumstances, this restriction is the true space of the problem and 
the solution is then exact. If we restrict |*,f > to the family of coherent statesfll, 12, 15-17] 
this corresponds, in ordinary quantum mechanics to a wave-packet of the form 



Here TV, A are the normalization and phase respectively and q, p specify a point in the classical 
phase space of dimension 2n. It is straightforward to show [18] that to 0(%), the equations 
for this system are Hamilton's equations for p, q, conserve the norm TV and yield 



where H is the classical Hamiltonian. Hence, \{t) is equal to the classical action; this can be 
used [18] to identify the semiclassical Bohr-Sommerfeld phase with the Aharanov-Anandan 
form of Berry's phase[19]. A more interesting situation arises when, as below, we consider 
squeezed coherent states (the wave-packets are allowed to spread). Along with the dynamics 
for the centroid of the wave-packet, we now also consider equations of motion for its spread. 
The full equations give us the mean- field theories mentioned above [9-12,20]. If, alternatively, 
in the Taylor series expansions around the centroid, we truncate to second order in deriva- 
tives of the Hamiltonian, we recover the equations of Heller [21], who has used this method 
extensively in studying 'quantum chaos'. 




(1) 




(2) 



A = p.q - iJ(q,p), 



(3) 
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Implicit in these 'semi- quantum' methods is the assumption that these dynamics are 
qualitatively similar to those of either the fully classical limit or the full quantum system. 
If the classical behavior is integrable, the semiquantal dynamics are said not to break the 
integrability[22]. Consideration of quantal fluctuations is supposed to 'supress' chaos and the 
full quantum dynamics are usually said to be regular. In this Letter, we demonstrate that 
this is not always true. Using the example of the double-well potential in ordinary quantum 
mechanics, we show that quantal effects, in fact, induce chaos. While neither the quantum nor 
the classical system are understood to display the sensitive dependence to initial conditions 
characteristic of chaos, the semiquantal dynamics does do so. We justify here the use of the 
terminology 'semiquantal'. This is not a trivial point - we derive the dynamics directly from 
the quantum Hamiltonian with no reference to the classical limit. In fact, it can be shown [22] 
that semiquantal dynamics exist even for systems without a well-defined classical dynamics. 
However, if the traditional 'semiclassical' and classical limits do exist, they can be recovered 
as special limits of semiquantal dynamics. The terminology is hence most apposite. 

While the standard approach uses the TDVP, we provide here an alternative derivation 
of semiquantal dynamics via the Ehrenfest theorem. We believe that this is somewhat 
more intuitive and is similar to formulations in non-equilibrium statistical mechanics. The 
dynamics obtained are the same as in the TDVP approach. Consider a particle of unit 
mass moving in a one-dimensional time-independent bounded potential with a Hamiltonian 
H = ^ + V(x ) where 6 denotes operators. The equations of motion for the centroid of a 
wavepacket representing the particle are 



where the <> indicate expectation values. In general the centroid does not follow the classical 
trajectory. We now expand the equations around the centroid using the identity 



where = d n F/du n \ <il> and U = u- < u > (the repeated index summation convention is 
used throughout, unless otherwise specified). Using this and operator commutation rules, 
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> 



(4a) 



(46) 



<F(u) >= < U n > F( n \ 
n 



n > 



(5) 



we can generate a countably infinite number of moment equations (corresponding to the 
infinite dimensional Hilbert space of the problem). The assumption that the wavepacket 
is a squeezed coherent state renders the space finite; it provides the relations < X 2m >= 
(2 ^^T (no summation), < X 2m+1 >= 0,4/i < P 2 >= h 2 + a 2 and< XP + PX >= a, which are 
easily recognised as those for generalized Gaussian wave-functions [9- 12, 20]. (This assumption 
is precisely that of the TDVP: the wavepacket is restricted to a given subspace). This yields 
the equations 

dx (a \ 

* = _jJLy(2™+i) (a; ) m = ,l,... (66) 
dt m\2 m w ' ' ' v ; 

- = a, (6c) 

g = g±^- 7 Q - 2 V^\x). m=l,2... (6d) 
dt 2p (m- l)!2 m - 2 w v ; 

The system is now reduced to the dynamics of x,p, p and a. ( where we write x,p for < x >, < 

j3 >) and are exactly those derived from the action principle [23]. 

We now introduce the change of variables p = p 2 , and a = 2pn, transforming equations 

(6) to 

dx \ 

^ (7a) 

dp = y(2ro+1) Q1 {n) 

dt m\2 m K h y 1 

<*> 

W-tt fi 2 n 2m-l 

^ = j-- - ^ ^(x), m=l,2... (7d) 

V (m- l)!2 m ~ 1 v ; v ; 

Remarkably, these new variables form an explicit canonically conjugate set, yielding a classical 

Hamiltonian phase space as our approximation to the Hilbert space. The classical degrees of 

freedom are the 'average' variables x,p and the 'fluctuation' variables p,n, respectively; the 

associated Hamiltonian is 

p 2 TT 2 

H ex t = y + y + V ext {x, p); (8a) 
fi 2 n 2m 

Vext ^p) = V(x) + — + J^V^Hx), m=l,2... (86) 

where the subscript ext indicates the 'extended' potential and Hamiltonian. This formulation 
is very interesting and possibly quite powerful; it provides us with an explicit gradient system, 
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and the extended potential provides a simple visualization of the geometry of the semiquantal 
space. We may thus get a qualitative feel for the semiquantal dynamics before proceeding 
to detailed (numerical) analysis. We note here that a) both the fluctuation and average 
variables are treated on the same footing and the phase space is dimensionally consistent: 
p has the dimensions of length and it that of momentum, b) the value of H ext is < H > 
under this approximation, and is conserved, c) V ext has an infinite barrier at p = 0, hence 
'quantum fluctuations' can not be zero except in the limit K — » and d) V ext is symmetric 
in p corresponding to a choice of sign in p = y/ji = V< X 2 >; the infinite barrier renders the 
choice one of convenience and of no physical significance. 

We note also that this formulation is exact for the simple harmonic oscillator. For this 
system, the variables (x,p) decouple from (p,n). The average variables hence execute the 
usual classical harmonic motion. The fluctation variables are a bound one-degree-of-freedom 
problem, and in general execute oscillatory motion. It is possible to find a fixed point in 
these variables, however. This is the familiar example of a Gaussian wavepacket that executes 
harmonic motion with a fixed spread (see Ref. 21). The extended potential, in this case, 
provides us with the quantum corrections due the 'zero-point fluctuations'. 

We consider now the dynamics of the simple double-well system with the Hamiltonian 
H = - 2|i + |i 4 . The extended Hamiltonian is 



where we have set h = 1. If we examine V ext in the upper half of the x - p plane, it is 
straightforward to show that there is always a minimum for the potential at (x,p) = (0,p), 
where p is the largest positive real solution of ap 4 - 3bp 6 + 1/4 = 0. For £ = < 1, this is the 
only minimum; quantum effects are so large that the the well barrier is effectively absent. 
As £ increases, two minima appear in the half-plane, corresponding to the well minima of 
the original problem; two saddle-points also emerge, interestingly. Also, the (x,p) system 
as driven by the (p,ir) system is like a non-linearly driven Duffing oscillator [24] with back 
reaction. With these ingredients, it is not surprising that the equations 




(9) 



dx 



(10a) 
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— = ax — bx — 3bxp , 
dp 



(106) 



(10c) 



§ = i+^-3^)-3^ 




(lOd) 



(which are the just the explicit form of equation 7 with F(p, £) as above) exhibit chaos. We 
have numerically investigated these equations over a range of parameter values and initial 
conditions. We summarize here some of these results at the typical parameter values of a = 
1,6 = .01 i.e., where quantum effects are small but not negligible. At low energies, as is typical 
for Hamiltonian systems, there are only integrable orbits. As the energy increases, chaotic 
orbits also emerge; the two types are very near each other in choice of initial conditions. The 
details of the phase-space structure will not be discussed here. However, the existence of 
chaos in this system has been established using Poincare sections and Lyapunov exponents 
calculations [25]. In Figs. 1,2 we show a typical Poincare section; taken with p = 0,p > 0. 
We can see the features characteristic of Hamiltonian chaos: the typical stochastic web 
structure, and in Fig. 2, which is an enlargement of Fig.l, the familiar appearance of islands 
in this stochastic sea [26]. The associated largest Lyapunov exponent X max equals 0.125. The 
calculation of Lyapunov exponents is also typical, with a slow saturation, and a residual 
oscillation of about 5%. The exponents are symmetric around and hence sum to (to 
within 10~ 6 ). We note that for various orbits (as for this one) where the classical motion 
would have remained confined to one well, the fluctuation dimension provides a means for 
the system to 'tunnel' between the two wells. 

In the classical limit, this system has only two independent first-order differential equa- 
tions, insufficient for chaos. The quantum dynamics are said to have 'no chaos', based on 
the fact that Schrodinger's equation is a linear partial differential equation; the dynamics 
are said to be formally equivalent to a classical harmonic oscillator system in quadrature 
on a 2iV-dimensional classical phase space. It has been recently argued[27] that this is too 
simplistic: the constraint of square integrability and arbitrary ray rotation renders the space 
a compact complex projected space CP(N-l) which has a different topological structure and 
may display substantially different behavior. There is a conjecture[22], further, that in sys- 
tems with inherently infinite-dimensional Hilbert spaces, semiquantal effects enhance chaos 
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(see Ref. 22 for details). For this system, the dynamical effect of the fluctuation variables is 
to induce chaos. This may be regarded as support for the arguments of Ref. 22; if, in fact, 
the full quantum system has 'no chaos', this is anomalous behavior that exists only in the 
semiquantal limit. 

In summary, we have introduced above a way of viewing the time-dependent Hartree- 
Fock method and other mean-field theories - the extended potential approach. Using this 
formulation for the double-well system, we have demonstrated that quantum fluctuations 
may induce chaotic semiquantal behavior, in keeping with a recent conjecture. Irrespective, 
care is indicated [28]: while variational and other classical approximations made to quantum 
systems are quite powerful, their qualitative behavior may be anomalous and may not exist 
in the full quantum system nor its completely classical limit. 

AKP acknowledges the Robert A. Welch Foundation (Grant No. F-0365) for partial 
support. 
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Figure Captions 

1. Fig.l. Poincare sections in the x,p plane at an energy of H ext = -1.25. There are about 
30000 points shown. \ max = 0.125. 

2. Fig. 2. Enlargement of Fig.l. Notice the distinctive stability islands in the stochastic sea. 
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